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Abstract 

A discrete stochastic process involving random amplification with additive noise is 
studied analytically. If the non-negative random amplification factor b is such that < 
b 13 >= 1 where (3 is any positive non- integer, then the steady state probability density 
function for the process will have power law tails of the form p(x) ~ This is a 

generalization of recent results for < (3 < 2 obtained by Takayasu et al. in Phys. Rev. 
lett. 79, 966 (1997). It is shown that the power spectrum of the time series x becomes 
Lorentzian, even when 1 < f3 < 2, i.e., in case of divergent variance. 
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Power law behavior of distribution function is widely observed in nature 0. Recently, 
Takayasu et al. presented a new general mechanism leading to the power law distribution 
P]. They analyzed a discrete stochastic process which involves random amplification 
together with additive external noise. They clarified necessary and sufficient conditions 
to realize steady power law fluctuation with divergent variance using a discrete version of 
linear Langevin equation expressed as 

x(t + l) = b(t)x(t) + f(t), (1) 

where f(t) represents a random additive noise and b(t) is a non-negative stochastic coef- 
ficient. They derived the following time evolution equation for the characteristic function 
Z(p,t) which is the Fourier transform of the probability density p(x,t): 

Z(p,t+1) = W(b)Z(bp,t)db$(p), (2) 

J 

where W(b) is the probability density of b(t) and $(p) is the characteristic function for 
f(t). They showed that when (b 13 ^ = 1 holds for < (3 < 2, the second moment {x 2 (t)} 
diverges as t — > oo, but Eq. (Q) has a unique steady and stable solution : 

lim Z(p, t) = Z(p) = 1 - const x \pf + ■ ■ • , (3) 

which yields the power law tails in the steady probability density 

lim p(x, t) = p(x) ~ l/x^ 1 , (4) 

t — >oo 

or equivalently, the cumulative distribution 

P(> \x\) ~ l/x p . (5) 

They also made numerical simulations of Eq. (|I|) by employing a discrete exponential 
distribution for W(b), and showed that the theoretical estimate of the relation between 
(3 and the parameters specifying W(b) (Eq. (15) in 0) nicely fits with the simulation 
"even out of the range of applicability, (3 > 2" . They state that "The reason for this lucky 



coincidence is not clear", although they point out at the same time that the power law 
distribution tails are a generic property of Eq. ([I]) 0. 

In this Brief Report, the following two statements will be presented: 



(A) Takayasu et al.'s theory can be straightforwardly extended for (3 > 2: If (b^J = 1 
holds for a positive non-integer (3, then there exists a unique steady and stable solution 
of Eq.® 

n 

Z(p) = J2 A 2m (-l) 2 V m ~C\pf + 0(p 2 "+ 2 ), (6) 

m=0 

where In is the largest even number that is smaller than j3. This Z(p) leads to p(x) ~ 



(B) When (bPj = 1 for a non-integer (3 between 1 and 2, the power spectral density 
(PSD) of x(t) is Lorentzian increasing with the observation time T as 
2 xl (1/Tj) (b 2 f 

S ^ T) ~f \n(b>) + ^ f ° r T>>1 ' (7) 

where 

4 = (AO)) + J^J, (8) 

Tl = In (6 2 ) + ln[l/(6)]' (9) 

From the statement (A), "the coincidence" found in is naturally understandable. 
To prove (A), we assume the following form for Z(p): 

oo oo 

Z{P) = E ^ + \pf E C np n , «0 = 1, (10) 

n=0 n=0 

and substitute it into Eq. (Q) in the limit t — > oo. If $(p) is an even function (i.e., the 
distribution function of f(t) is symmetric as assumed in 0), we can first prove that a\ = 
because (b) ^ 1. Also, ci = because (fr^ 1 ) ^ 1. Thanks to a 2m -i = and (b 2m+1 ) ^ 1, 
o-2m+i = is derived. Similarly, c 2m _i = and /y 3 + 2m + 1 ^ ^ i yield c 2m +i = 0. We can 
thus prove that a n and c n in Eq. flTU| ) vanish for all odd numbers n, i.e., Eq. (|B|) holds. 
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(Note that the n th moment (x n (t)) with n > (3 diverges not only for even number n but 
also for odd number which corresponds to the vanishing coefficient a n .) Taking exactly 
the same procedures as in ||, we can prove that this solution is unique and stable. In 
case of (3 > 2, we have a finite variance but higher order moments, (x n (t)) with n > (3, 
diverge as t — > oo. 

To derive the probability density p(x), we only need to assume that all k-th derivatives 
of Z(p) satisfy the boundary condition 

lim d k Z{p)/dp k = 0. (11) 

p— >±oo 

Using Eq. ([O]), we can partially integrate the expression 

p{x) = ^J™J xp Z{p)dp (12) 

\P] + 1 times, where \0\ is the largest integer that is smaller than (3. Thus we obtain the 
asymptotic expansion as 

p(x) ~ / e-*\Zf-W- l dt ~ \x\- {(3+1) T{(3 - [(3}), (13) 



where T is the Gamma function. 

To prove the statement (B), we note that the two-time correlation function is rigorously 
obtained from Eq. (|T|): 

<f>(T,t) = (x(t + T)x(t)) = (x 2 (t))(b) T , (14) 

where 

(*t)) = (P)\*{0)) + \d&(f'). (15) 

If 1 < (3 < 2, we have a relation < (b) < 1 < (b 2 ) because the function G(j) = (b 1 ) 
satisfies G(0) = 1 and G"{^) > @. Then (ft increases with t, but decays with r as 
~ e~ T / T ° for any fixed value of t (Debye-type relaxation), with the relaxation time 



Since the correlation function depends on both r and t, the Wiener-Khinchin relation 
cannot be used to obtain the PSD. Defining the PSD which depends on the observation 
time T as 



and using 4>(r,t) obtained above, we arrive at the expression (|7|). The spectrum is 1/f 2 
-type for / ^> 1/t\ and flat for / <C 1/Vi. Equation ([!]) implies that the power increases 
exponentially with the observation time T, which corresponds to the divergent behavior 
of the variance (x 2 (t)). (We have neglected the case < (3 < 1 where even the average of 



When (3 > 2, both < (b) < 1 and < (b 2 ) < 1 hold and results are rather trivial: 



Thus, as far as the PSD is measured, we cannot observe any singular aspect, higher order 
singularities being hidden. 

The stochastic process described by Eq. (1) generally leads to the power law behavior 
p(x) ~ l/x@ +1 , while it also yields a Lorentzian spectrum S(u>) oc l/[(l/r) 2 + u 2 ]. A 
colored noise, or l/f a fluctuation, whose PSD is proportional to l/uo a has attracted 
much attention since 1/f noise was discovered several decades ago||. Such power law 
behavior of the PSD is also observed widely in nature, and these two power laws, one in 
the probability density and the other in the PSD, are sometimes discussed together [2]. 
Therefore it is interesting to know whether an extremely long time scale r can be involved 
in the present stochastic process. Because, in that case, the observation time T, which 
relates with the low frequency cut-off ujq as ujq = 2n/T, cannot reach this time scale, then 
a l// 2 fluctuation, namely S(uj) ~ 1/u) 2 (for uj > uj ) is observed practically. 




(17) 



x diverges as (x(t)) = (b) (x(0)) because (b) > 1.) 




(18) 



(19) 



(20) 
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One can immediately see that the time constant t or T\ becomes large in very limited 
cases. First, the average of b should be close to unity, i.e., (b) = 1 — e with < £ « 1. 
Then To becomes ~ l/e> 1. Furthermore, in case of (5 > 2, we need (b 2 ) smaller than 
unity, while in case of 1 < f3 < 2, the condition (b 2 ) = 1 + 5 with < 5 <C 1 is necessary. In 
the latter case, we obtain T\ ~ l/(e + 8) ^> 1. The exponential or Poisson distribution for 
W(b) does not lead to such a long time constant. One example of large T\ is obtained by 
choosing W(b) to be a narrowly peaked distribution having the average which is slightly 
smaller than unity and the second moment slightly larger than unity. 

As pointed out in the above, it should be noted that a stochastic process whose 
stationary density function has power law tails will not necessarily exhibit the power law 
behavior in PSD. 
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